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1 Introduction and summary 



Pion-nucleon scattering is an important testing ground for our understanding of the chiral dynamics 
of QuantumChromoDynamics (QCD). It has been investigated to third and fourth order in the chiral 
expansion, leading e.g. to precise predictions for the threshold parameters [1, 2]. This approach 
is based on an effective field theory with the active degrees of freedom being the asymptotically 
observable pion and nucleon fields. When going to higher energies, the usefulness of the chiral 
expansion is limited by the appearance of the nucleon resonances, the most prominent and important 
of these being the A(1232) with spin and isospin 3/2. Its implications for hadronic and nuclear physics 
are well established. Consequently, one would like to have a consistent and systematic framework 
to include this important degree of freedom in baryon chiral perturbation theory, as first stressed 
by Jenkins and Manohar [3] and only recently formalized by Hemmert, Holstein and Kambor [4]. 
Counting the nucleon-delta mass splitting as an additional small parameter, one arrives at the 
so-called small scale expansion (which differs from the chiral expansion because the NA splitting 
does not vanish in the chiral limit). It has already been established that most of the low-energy 
constants appearing in the effective chiral pion-nucleon Lagrangian are saturated by the delta [5] 
and thus the resummation of such terms underlying the small scale expansion (SSE) lets one expect 
a better convergence as compared to the chiral expansion. In addition, the radius of convergence is 
clearly enlarged when including the delta as an explicit degree of freedom. This of course increases 
the complexity of the approach since the pertinent effective Lagrangian contains more structures 
consistent with all symmetries. The purpose of this paper is to analyze pion-nucleon scattering to 
third order in this framework. In particular, we want to address two questions. First, it has to be 
demonstrated that for a given energy range, the third order SSE calculation leads to improved results 
as compared to the ones obtained in the third order chiral expansion. Second, a precise description 
of the resonant phase, i.e. of the P33 partial wave, should be obtained. As we will demonstrate, the 
SSE will pass both these tests. Furthermore, our investigation can be used to test the convergence of 
the small scale expansion. Indeed, from many models and more phenomenological approaches it is 
believed that the most important delta contributions stem from the Born graphs with intermediate 
resonance fields. We will address this issue in what follows. Also, from the technical point of view, we 
provide for the first time a systematic evaluation of the many 1/m corrections (using the machinery 
spelled out in Ref.[4]). Finally, we wish to point out that a systematic Lorentz invariant formulation 
as it is available for the pion-nucleon effective field theory [6] is not yet available, but could be built 
based on the pioneering investigation in Ref . [7] . The heavy fermion approach used here suffers from 
the same deficiencies as the standard heavy baryon chiral perturbation theory, these problems are, 
however, not relevant to the main issue we are going to address, namely the extension to higher 
energies. 

The pertinent results of the present investigation can be summarized as follows: 

(i) We have constructed the one-loop amplitude for elastic pion-nucleon scattering based on an 
effective field theory including pions, nucleons and deltas to third order in the small scale 
expansion, 0(e 3 ), where e collects all small parameters (external momenta, the pion mass and 
the nucleon-delta mass splitting). We have constructed the pertinent terms of the effective 
Lagrangian including the 1/m corrections. The amplitude contains altogether 14 low-energy 
constants from the nucleon, the nucleon-delta and the delta sector (if we count the leading 
7riVA coupling constant g n NA as a LEC). 
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(ii) The values of the LECs can be determined by fitting to the two S- and four P-wave amplitudes 
for different sets of available pion-nucleon phase shifts in the physical region at low energies. 
We have performed two types of fits. In the first one, we fit to the Matsinos phase shifts in 
the range of 40 to 100 MeV pion momentum in the laboratory frame. This allows for a direct 
comparison with the results based on the chiral expansion. We find that the third order SSE 
results are clearly better than the ones of the third order chiral expansion and only slightly 
worse than the ones obtained at fourth order in the pion-nucleon EFT. Second, we have fitted 
to the Karlsruhe phase shifts for pion lab momenta between 40 and 200 MeV. This allows for 
a better study of the resonance region. In both cases, most fitted LECs are of "natural" size. 

(iii) We have studied the convergence of the small scale expansion by comparing the best fits based 
on the first (leaving the coupling constant g n NA free), second and third order representation 
of the scattering amplitudes. The third order corrections are in general not large, but they 
improve the description of most partial waves. This indicates convergence of the small scale 
expansion for this process. As anticipated, the most important contributions come from the tree 
(Born) graphs with intermediate delta states. This allows to pin down the coupling constant 
g wNA fairly precisely. 

(iv) We can predict the phases at /operand at higher energies, in particular the threshold parameters 
(scattering lengths and effective ranges). The results are not very different from the third and 
fourth order studies based on the chiral expansion, but the description of the scattering length 
and the energy dependence in the delta channel are clearly improved. While the convergence of 
the isovector S-wave scattering length is satisfactory to third order in the small scale expansion, 
for drawing a conclusion on the isoscalar S-wave scattering length a fourth order calculation is 
mandatory. 

(iv) We have considered the pion-nucleon sigma term. To third order in the small scale expansion, 
it depends on the low energy constant c\ and the coupling g n NA- For the KA85 phases, we get 
a value consistent with previous determinations. We have also performed fits using the sigma 
term extracted from a family of sum rules as input and shown that this procedure leads to 
more reliable predictions for some LECs. 

The manuscript is organized as follows. In section 2, we discuss the effective Lagrangian underlying 
our calculation. It decomposes into separate nucleon, delta-nucleon and delta contributions. We 
spell out all details pertinent to our calculation including also the 1/m corrections. Many of these 
terms have so far not been available in the literature. Section 3 contains the results for the pion- 
nucleon scattering amplitudes g ± ,h ± to third order in the small scale expansion. In particular, we 
discuss the interplay between the low-energy constants from the various sectors and the role of the 
so-called off-shell parameters. The fitting procedure together with the results for the phase shifts 
and threshold parameters are presented in section 4. In particular, we give a detailed comparison 
with the results obtained in the chiral expansion and discuss the issue of convergence. We also discuss 
the sigma term and show how it can be used to further constrain the fits. The pertinent Feynman 
rules, explicit expressions for the various 1/m corrections, for the tree and loop contributions to the 
scattering amplitude and the analytical expressions for the various threshold parameters are given 
in the appendix. 
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2 Effective Lagrangian 



In this section, we briefly discuss the effective Lagrangian underlying our calculation. We borrow 
heavily from the work of Hemmert et al. [4] and refer the interested reader to that paper to fill in the 
details omitted here. The explicit inclusion of the delta into an effective pion-nucleon field theory 
is motivated by the fact that in certain observables this resonance plays a prominent role already at 
low energies. The reason for this is twofold. First, the delta-nucleon mass splitting #4 is small, 

A = m A ~m N = 294 MeV ~ 3F n , (2.1) 

with F n = 92.4 MeV the weak pion decay constant. In the chiral limit of vanishing quark masses, 
neither A nor F n vanish. Therefore, such an extended EFT does not have the same chiral limit as 
QCD, as it is well-known since long [8]. Second, the delta couples very strongly to the irN^ system, 
e.g. the strong AV7 Ml transition plays a prominent role in charged pion photoproduction. One 
can set up a consistent power counting by using the well-known heavy baryon techniques [9, 10] 
and by treating the mass splitting A as an additional small parameter besides the external momenta 
and quark (meson) masses. Therefore, any matrix element or transition current has a low energy 
expansion of the form 

M = e n Mi + e n+1 M 2 + e n+2 M 3 + 0(e n+3 ) , (2.2) 

where the power n depends on the process under consideration (for pion-nucleon scattering, n equals 
one) and e collects the three different small parameters, 



^ f.^r . ( 2 - 3 ) 




with A x ~ 1 GeV the scale of chiral symmetry breaking, M n the pion mass, and q some external 
momentum. This power counting scheme is often called e- or small scale expansion (SSE). As it 
is common in heavy baryon approaches, the expansion in the inverse of the heavy mass and with 
respect to the chiral symmetry breaking scale are treated simultaneously (because rriN ~ rn& ~ A x ). 
The effective Lagrangian has the following low-energy expansion 

£ eff = + £ {2) + £ (3) + . . . (2.4) 

where each of the terms decomposes into a pure nucleon (irN), a nucleon-delta (irNA) and a 
pure delta part (it A), 

£W = £S + + 42 , n = 1, 2, 3, . . . . (2.5) 

We will now discuss these in succession. We note that the NA terms have to be understood symmet- 
rically, i.e. we can have an in-going nucleon and an out-going delta or the other way around. The 
coupling of external fields is done by standard methods, for our purpose we only have to consider a 
scalar source to deal with the explicit chiral symmetry breaking due to the quark masses. 

# 4 We use the same symbol for the delta resonance as well as the AN mass splitting. This cannot lead to confusion 
since from the context it is always obvious what is meant. 
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2.1 Single nucleon sector 

This part is fairly standard, for completeness we collect the terms necessary for the following discus- 
sion. To lowest (first) order, the relativistic pion-nucleon Lagrangian takes the form 



4£ = 



(2.6) 



where the bi-spinor \I/ jv collects the proton and neutron fields, qa = 1-26 is the axial-vector coupling 
constant and m = itln the nucleon mass. All parameters in the Lagrangian should be taken at 
their chiral limit values. We do not exhibit this by special symbols but it should be kept in mind. 
The heavy baryon projection is most economically done using path integral methods as outlined 
in Ref.[10]. This leads to a set of matrices, which organize the transitions between the light-light 
(N — N), light-heavy (N — h) and heavy-heavy (h — h) components of the nucleon fields, denoted 
by A, B and C, in order, 



C^n = NA N N + hB N N + N lo B^ Q h - hC N h 



(2.7) 



The inverse of C is then further expanded in inverse powers of the nucleon mass, leading to the 
decoupling of the positive and negative velocity sectors. These matrices have a chiral expansion, i.e. 
A = A^ + A^ + . . . . From the lowest order, we only need 



A {1) 



B 



(i) 

N 



Ai) 



2m , 

iv ■ D + g^S ■ u , 
- l5 (2tS-D + ^vu) , 

iv ■ D + qaS ■ u , 



(2.8) 
(2.9) 
(2.10) 

(2.11) 



in terms of the nucleon four-velocity t> M and the Pauli-Lubanski spin-vector (for more details, 
see e.g. the review [11]). We now turn to the second order terms. In the isospin limit of equal up 
and down quark masses, one has to deal with 4 operators, 



,(2) 

'ttJV 



N 



ci(x+) ~ ^({u,u u ){D^D u } + h.c.) + |( M 2 ) + 



c 3 



?C4 



N 



(2.12) 



where x+ includes the explicit chiral symmetry breaking and traces in flavor space are denoted by 
(. . .). More precisely, x± = u^X^ =t u X^ u i u n — i^d^u — ud^v)) and U(x) = u 2 (x) collects the 
pion fields (for details, see e.g. [11]). Note, however, that due to the presence of the delta degrees of 
freedom the numerical values of these LECs are different from the ones obtained in the pure pion- 
nucleon EFT. This is discussed in more detail in Ref.[12]. From the second order transition matrices 
we need the following terms: 



A {2) 

»(2) 



ci(x+) + c 2 (v ■ uf + c 3 u 2 + c 4 [5' /1 , S^u^Uy , 
-c 4 7 5 [w • u, S ■ u] , 

- ci(x+) + c 2 {v ■ u) 2 + c 3 u 2 + c 4 [S ,/1 , S u ]u^u v 



(2.13) 
(2.14) 

(2.15) 



4 



Finally, from the third order Lagrangian, we only need one 1/m correction to the dimension two 
operator ~ c 2 and the following dimension three operators 

A " = 7^(( v - u K) DI * + h - c -) +idiK,[v-D,u^]+id 2 [u^[D^vu}} 
+id 3 [v ■ u, [v • D, v • u}} + d 5 [x-,v ■ u] - i [S^, S u ]d u ([v • D,u^}u u ) 
-i \S», S u ]d lb {u^[D u , vu}) + d 16 ( X+ )S -u + td 18 [S-D, x-] ■ (2.16) 

Note that in nN scattering, some of the LECs only appear in certain combinations, here d 1 + d 2 and 
di4 — d±5 are of relevance. Some of these LECs are needed for the renormalization and their finite 
parts depend on the regularization scale. We do not further specify this but it should be kept in 
mind. Furthermore, there are some additional LECs just needed for the renormalization, as spelled 
out e.g. in Ref.[l]. We refrain from writing down such terms here. As stressed before, the numerical 
values of the finite parts of these LECs are influenced by the presence of the delta and thus can not 
be taken over from the pure 7rN EFT. The remaining 1/m corrections from the elimination of the 
small components are standard and can be found in Ref. [1]. 



2.2 NA-sector 



We now turn to the sector of the nucleons coupled to deltas and pions as well as external sources. 
Here, we only consider external scalar sources related to the explicit chiral symmetry breaking. The 
first order relativistic Lagrangian consistent with the requirement of point transformation invariance 
takes the form 



c 



with & ail (z 



W - n 



(2.17) 



9ix,v + ZqIhIv, ^ is a conventional Rarita-Schwinger spinor and w l a = j(r % u a ). 
We have to deal with two parameters, the leading pion-nucleon-delta coupling constant g n NA and 
the so-called off-shell parameter zq. The latter will be discussed in more detail below. The heavy 
baryon projection is more tedious since the delta field has a large and a small spin-3/2 as well as 
four (off-shell) spin-1/2 components. To keep only the large (light) spin-3/2 part, one has to insert 
appropriate spin-isospin projection operators. The technology to do that is spelled out in Ref. [4]. 
The effective Lagrangian has the genuine form 



£nNA = TA NA N + GB NA N + T^V ] N/ ^ h + G^C ] NA ^h + h.c. 



(2.18) 



i.e. one has to deal with four types of transition matrices. These are denoted by Ana (light 
nucleon (N) to light-delta transitions (T)), Bna (light-nucleon (N) to heavy-delta transitions (G)), 
Cna (heavy-delta transitions (G) to heavy-nucleon (h)), and V NA (light-delta transitions (T) to 
heavy-nucleon (h)). Note that the field G has 5 components. To order e, these transition matrices 
read 



NA 



B 



(i) 

JVA 







9wNA 



-*&f&P+S ll S-w i P+ 
2z P^ 5 v^S ■ w l P + 
-2z Q P- lb S lx vw i P + 
\ (1 + Z )P+Vf,V ■ w i P + 



(2.19) 



(2.20) 



5 



V 



(1) 

NA 



o, 



^JVA ~~ i/7riVA 



P_< 3 P a/1 P- 

2z P_t; • 11^75 P+ 
(1 + z )P_?; • w'vpP- 

-4(1+3*,,) p_ i g . 

V -2z P_S • w^7 5 P+ / 



r2.211 



(2.22) 



with 



4: 



(2.23) 



in d space-time dimensions and the P± are the usual velocity projection operators. The second order 
relativistic nNA Lagrangian has two terms of relevance to our study, these read 



+ i°A w i RiD e 

m 



^iv + h.c. , 



(2.24) 



with z another off-shell parameter also discussed below. To be more precise, any new structure 
which appears with a low energy constant has a separate off-shell parameter. However, these can 
be absorbed in the corresponding LECs (as discussed in more detail below) and we thus collectively 
call these new off-shell parameters z. The LECs hi (i = 3, 8) are finite and will appear in the tree 
contribution to the ttN scattering amplitude. The heavy baryon projection is done in terms of the 
appropriate second order transition matrices 



A {2) 



B 



NA 

(2) 

JVA.l 



P+ 3P^J(h + b 8 )w t a/3 v l3 , 
-P_ 3 P^ 5 2ib 3 w i a(3 S p P + 



B 



(2) 

NA,2 



a ps 

3 P+^ 



B 



(2) 

JVA,3 



B 



(2) 

NAA 



B 



(2) 

ATA, 5 



V 



(2) 
NA 



M3) 



1 + 3z)i(b 3 + b a )w % a() S a v p - 3zib 3 w l ap v a S p 
P-^75 [ - (1 + z)2ib 3 w i ap v a S^ + 2zi(b 3 + b 8 )w i a/3 S a v p 

4 r -3 

-P_S , M 7 5 [(1 + 3z)2ib 3 w l al3 S a S l 3 - -zi(b 3 + b 8 )w t a/3 v a v^ 
P+Vp[(l + z)i(b 3 + b 8 )w l al3 v a vp - 2zib 3 w l al3 S a S p 

P-l&ibzSffljpa 3 P a ^P+ , 

ib« 



NA 



P+ ?,PfJLa 



m 



w l a/3 iDp . 



(2.25) 
(2.26) 

(2.27) 
(2.28) 
(2.29) 
(2.30) 
(2.31) 
(2.32) 



Note that at this order the first non-vanishing contribution to T> N & appears. Finally, at third order, 
we only need the relativistic effective Lagrangian and the corresponding transition matrix A$a- 
These read: 



£Sa =%®»>M 



^iv + h.c. , 



(2.33) 
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= i-w)JD v + (/i + / 2 )[L> M ,< /3 ] W + / 4 <(x + > + h\D^X 1 -} • (2.34) 

Four new LECs, which we call fi, appear. However, only the combinations f\ + fi and 2/ 4 — / 5 are of 
relevance in case of pion-nucleon scattering. Altogether, in the NA sector we have 5 LECs and the 
same number of (unobservable) off-shell parameters. We have not counted 63 and bs separately, since 
the latter can be absorbed in other LECs as detailed below. The Feynman rules for the resulting 
7riVA vertex to first, second and third order are collected in appendix A. The 1/m corrections 
originating from the elimination of the various "small" components are discussed below. 



2.3 Single A sector 

As in the previous paragraphs, we start with the dimension one effective Lagrangian coupling the 
massive spin-3/2 fields to pions, 

+f g m ^75 + |(7X j + u% lv )^ + | 7 m fit-mM , (2-35) 

with u l i = 5^ and D l i the appropriate chiral covariant derivative. If one now performs simultane- 
ously the 1/m expansion of the nucleon and delta fields, one can only rotate away the nucleon mass, 
so that the NA mass splitting remains in the delta propagator. Since this is, however, a quantity 
of order e, it can be expanded systematically. In the heavy fermion approach, the Lagrangian takes 
the form 

ArA = TA A T + GB A T + T-yoBl-yoG - GC A G . (2.36) 

The corresponding contributions from the three types of transition matrices A A , B A , and C A are 
given by 

A W = -( iv . - AS^ + 9l S ■ , (2.37) 

Bfgi = P- (2iS ■ £>« + • ««) 76 P + , B« = -P + + g 2 )S^ P + 

= P - f V4?75 P + , B$ A = ^P^iD^ P + , B% = , (2.38) 
/»-i(o) _ 1 P P P r -i(o)_ 1 p p p r -i(o)_ 1 2 

L A,11 - "^m^- 3^^- ' L A,22 - 2^"^+ l^-M" ' L A,23 - _ 2m3 ^ ' 



^A,32° — ^P- v iiSvl5P+ j C A ^3 — ~1^P- 2P^yP~ , C A ]te ~ ~ 2m^~ iP^yP- , (2.39) 



= (^J P- 3 P^(^-^ + A^+ 5l( 5-^) 3 P Qi ,P_, 

Ca'i? = (^) 2P_ 3 P lia iS.D">y 51 P oa ,P + , 
-i(i) / 1 \ 2 2 



1 



2 
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/»-!(!) 
U A,14 

1(1) 



c 



A, 15 



^-1(1) 
"A, 22 



/»-!(!) 
U A,23 

/»-!(!) 
U A,24 



/»-!(!) 
U A,25 

/»-!(!) 
U A,33 



c 



-1(1) 

A, 34 



/»-!(!) 
U A,35 

/»-!(!) 
U A,44 



c 



c 



-1(1) 

A,45 

-1(1) 
A, 55 



( 1 \ 2 2 

< 1 (91 + ^92)P- sP^ a U^ 5 V u P + , 



\2mJ 
\2m) 3 

2 



-P+S^-iv ■ D« + 2AS" - -S ■ u«( 9l - 4 92 + 8g 3 )]S„P + , 



(A 

\2m/ 



2 • D" - \v ■ u ij ( 9l + 2g 2 )] l5 S v P- , 
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1 



\2mJ 9 
1 ^ 2 



-P+S^mS -D t3 +v u t3 (U 9l + 16g 2 + 16g 3 )]v u P + , 



V2m/ 

(w 
(2m) 



2 4 ID 

P-v^-iv ■ D ij + -AS ij + — 9l S ■ u l3 ]v u P_ , 

2 40 2 

1 



P-V^AiS ■ D" + -v ■ U v( 9l - 2g 2 )] l5 v u P + , 



(2m) 3 



" *P-S„[iv ■ D ij - 2A5 ij + \ 9l S ■ u ij ]S v P- , 



{.2m) 
(2m) 



2 4 1 fi 8 

P_^[--^ • £>« + -A6 i3 + S ■ u i3 \-A 9l + -g 2 )h 5 v u P + , 

2 P+v^-^iv ■ D ij + AA5 ij + ■ u ij (-17 9l + 12g 2 - 8g 3 )]v u P+ 



(2.40) 



The other elements of C A X are given by 



Ca^ = 7oCa;J7o • (2.41) 

From the second order relativistic Lagrangian, we only need one term, which is the mass insertion 
on the A propagator, i.e. the analog to the dimension two nucleon term ~ ci, 



41 = %e»Azi)*i(x + )s t3 9^®v>M)vl' 



(2.42) 



which translates into 



B 



(2) 



P + 3 P^a 1 (x+) g^'v'8 13 3 Pu'uP+ , 
. 



(2.43) 
(2.44) 



For our purpose, no dimension three operator of the nA sector is needed, thus 



(2.45) 



The genuine 1/m corrections will be considered in the next section. 



8 



2.4 1/m corrections 



To calculate the 1/m corrections related to the elimination of the "small" components, it is most 
economical to use the path integral formalism outlined in Ref. [10]. Here, an algebraic complication 
due to the multiplication of degrees of freedom to be eliminated and also due to the transition between 
the nucleon and the delta sectors appears. The machinery to do that is spelled out in great detail 
in Ref. [4], here we only collect the pertinent results and basic definitions. We note however that we 
give the most detailed list of such terms so far worked out, only a subset of these is found in Ref. [4]. 
After the canonical change of variables to make the action quadratic in the fields, the various pieces 
of the effective Lagrangian take the form: 



£>nN 
£wNA 

employing the definitions: 



N 



NA N N + N[i Q & Nl{ £^B N + i B na 1vC a 1 B N a 
TA a T + T[ 1o B a1o C a 1 B a + 1qV na1q C n 1 V N a]t , 
TA NA N + f\ lQ V NAlo C N 1 B N + toBItoCZ^jva! N + h.c. , 



B N 



Bn + CnaP a Bna , 
C-n ~~ C-naC- a C na i 
£Va + C N aC a x B a . 



(2.46) 
(2.47) 
(2.48) 



(2.49) 
(2.50) 
(2.51) 



We note that the only 1/m corrections for C kA we need are propagator insertions. From these 
equations our previously made remark that the low energy constants of the single nucleon sector 
get modified due to the presence of the delta becomes quite obvious since the light nucleon to light 
nucleon Lagrangian is modified by the appearance of the nucleon to delta transition operators, e.g. 
the last term in Eq.(2.46). We come back to this point below. The algebra to work out these terms 
is somewhat tedious and we collect the final results in appendix B. 



3 Pion— nucleon scattering 
3.1 Basic definitions 

In this section, we only give a few basic definitions pertinent to elastic pion-nucleon scattering. For 
a more detailed discussion, we refer to Ref.[l]. In the center-of-mass system (cms), the amplitude 
for the process 7r a {qi) + N(pi) — > it b {q 2 ) + N(p 2 ) takes the following form (in the isospin basis, with 
a, b denoting cartesian isospin indices of the pions): 



ivN 



E + m 
2m 



+i e bac r c 



g + (u, t)+ia- (q 2 x q\ ) h + (u, t) 
g~{u;,t) +ia- (q 2 x q 1 )h~(u;,t) 



(3.1) 



with uj 
and q \ 



v ■ qi = v 

q 2 = i = 



q 2 the pion cms energy, E\ = E 2 
[{s -Ml- m 2 ) 2 - Am 2 Ml] /(As), t 



E = (q 2 + m 2 ) 1 / 2 the nucleon energy 
- (qi — q 2 ) 2 is the invariant momentum 



transfer squared and s denotes the total cms energy squared. Furthermore, g ± (uj,t) refers to the 
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isoscalar/isovector non-spin-flip amplitude and h^iuj, t) to the isoscalar/isovector spin-flip amplitude. 
This form is most suitable for a heavy fermion calculation (as done here) since it is already defined 
in a two-component framework. 

The quantities of interest are the partial wave amplitudes /i±(s), where / refers to the orbital angular 
momentum, the superscript '±' to the isospin (even or odd) and the subscript '±' to the total angular 
momentum (j = I ± s), are given in terms of the invariant amplitudes via 



E + m f+ l 



dz 



g ± P l {z) + fh ± {P l±1 {z)-zP l {z)) 



(3.2) 



where z = cos(#) is related to the scattering angle. The Pi(z) are the conventional Legendre poly- 
noms. For a given isospin /, the phase shifts 6f±(s) can be extracted from the partial waves via 



2i\q\ 



exp(2^ 7 ± ( S )) - 1 



(3.3) 



For vanishing inelasticity, which is the case for the energy range considered in this work (y/s < 
1.3 GeV), the phase shifts are real. They are given by 



arctan(|g | Re fj±(s)) . 



(3.4) 



In the low energy region, one could equally well use the definition without the arctan, the difference 
being of higher order. For the phase shifts in the kinematical region considered here, this difference 
is not negligible. In fact, this arctan prescription is nothing but a unitarization procedure which is 
mandated by the appearance of the poles in the delta propagator at uj = A. Although there is some 
arbitrariness in this unitarization procedure, it has already been shown in Ref. [13] that it leads to 
the proper resonance width (for sharp resonances like the delta or the p in pion-pion scattering, for 
more details see Ref. [13]). Consequently, the phase shifts presented in what follows are based on 
Eq.(3.4). 



3.2 Small scale expansion of the amplitudes 

In this section we discuss the small scale expansion of the non-spin-flip and spin-flip amplitudes 
g ± ) h ± . These consist of essentially three pieces, which are the tree contributions with intermediate 
nucleons and deltas, counterterm parts of polynomial type as well as the unitarity corrections due 
to the pion loops (again with intermediate nucleons and deltas). The tree plus counterterm and 
loop graphs with intermediate deltas are shown in Fig.l and Fig. 2, respectively. We remark that 
the mrNA vertices only start at order 0(e 2 ) and thus all loop graphs of relevance here have NA 
couplings linear in the pion field. Formally, the small scale expansion of the various amplitudes has 
the form 

X = X trcc + X ct + X loop , X = g ± ,h ± , (3.5) 

where the tree contribution subsumes all Born terms with fixed coefficients, and the counterterm 
amplitude the ones proportional to the dimension two and three LECs. The last term in Eq.(3.5) 
is the leading one-loop amplitude consisting of terms of order e 3 . The latter is a complex- valued 
function and restores unitarity in the perturbative sense. Its various terms are all proportional to 
Note that we treat the chiral symmetry breaking scale A x ~ 1 GeV on the same footing as 
the nucleon and delta mass. These amplitudes are functions of two kinematical variables, which we 
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choose to be the pion energy and the invariant momentum transfer squared, i.e. X = X(u,t). In 
what follows, we mostly suppress these arguments. It is also instructive to compare these amplitudes 
with the ones of the pure pion-nucleon EFT, based on the chiral expansion. These terms are, of 
course, contained in the SSE to third order and their explicit expressions are given in Ref. [1] (we 
remind the reader that the numerical values of the LECs are different, for the reasons discussed 
above).* 5 The novel tree and loop terms, shown in Fig.l and Fig. 2, are of fourth order (or higher) in 
the chiral expansion, since the delta is frozen out in that approach and generates second (and higher) 
order contact interactions (as discussed in more detail in Ref. [5]). Therefore, the amplitudes based 
on the third order small scale expansion can be considered as partial fourth order chiral calculations. 
Of course, they also contain higher order terms in the chiral expansion, since the delta is not frozen 
out. One can therefore expect that in most channels the third order SSE calculation should lead 
to a better description of the phase shifts than given by the third order chiral expansion. Such an 
expectation is indeed borne out by the actual results to be presented below. 

The full one-loop amplitude to order e 3 is obtained after mass and coupling constant renormalization, 



(9a, rh, A, F, M) -> (g A , to, A, F n , M n ) . 



(3.6) 



The pion mass, decay constant and Z-factor are standard, 



o A 2M 2 A 

Ml = M 2 h+—£ 3 + 



2F 2 



1 _ 2M 2 



(3.7) 

(3.8) 
(3.9) 



in terms of the LECs £ 3i4 and the divergent pion self-energy tadpole A^ given e.g. in Ref. [11]. 
Similarly, for the nucleon mass shift, the nucleon Z-factor and the pion-nucleon coupling, we use 
(for more details, see Refs.[2, 12]): 



TO 



J N 



9A 



to -4M 2 ci + 



9 (9a\ 



4 \F J 



J 2 (0) 



+ 



4 (g 



tNA 



2 r 



(M 2 -A 2 )(J (-A) 
9 (9a^ 2 



2A 

M 2 

2 



A + l8^ (3M2 - 2A2) 



(3.10) 



4 / ^jva V 
3 V F J 



(M 2 -A 2 )J^(-A)+2AJ (-A) 
+ A(M 2 - 3A 2 )A. + — (— M 2 + 2A 2 ) 



M 2 



16tt 2 
K ~ 4n 2 



(3.11) 



# 5 In this paper, we only display the novel delta contributions to the amplitudes, threshold parameters and so on. 
The purely nucleonic terms can be found in Ref. [1] . 
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+ 



32 



I ( 4 _ 100 ^ 



81 9a, 



(M - A ) Jq(-A) + 2AJ (-A) 



27A 
1 



MVo(O) - (M - A ) J (-A) 



+ 



27M 2 ^ 



(M 2 -19A 2 )--^(M 2 -3A 2 



, V -M 2 + 2A 2 )(2 - — ^-) + -A-(-3M 2 + 2A 2 ) 



(3.12) 



in terms of the standard loop functions Jo^ij^) [H]- Finally, for the delta, we only need the mass 
shift, 

m A =m A -4MV - ^ (|) 2 M 2 J o (0) - ± (^)' (M 2 - A 2 ) (j (A) + ^A.) . (3.13) 

We refrain from giving the much more complicated expressions for the delta Z-factor and the renor- 
malization of g^NA- Also, we always work with renormalized LECs, all infinities appearing in the 
loop diagrams are accounted for by the corresponding infinite parts of the pertinent counterterms. 
What is still missing in the SSE is a systematic investigation of renormalization as it is the case in 
(heavy) baryon chiral perturbation theory. 



3.3 Counterterm combinations and off— shell parameters 

First, we enumerate the novel low energy constants related to the 7rA and tcNA sectors. Consider 
first the tiNA couplings. At leading order, there is only g n NA — 1-05 from the width of the decay 
A — > Ntt [4] . Since the necessary resummation at the pole of the delta includes some higher order 
terms, we will also perform fits leaving this coupling free. We expect, however, that the so determined 
value is not very different from 1.05. In addition, from the dimension two and three Lagrangians, 
we have the LEC combinations 63 + bs, fi + j'2 and 2/4 — f$. We note that there is a 1/m correction 
to b 8 , nonetheless the LECs 63 and bs need not be treated separately because the term bs /m can 
be absorbed in the dimension three LECs from the nucleon sector, as detailed below. From the An 
sector, we only have the coupling gi, which we leave free. In the large N c limit of QCD, one obtains 
the relation g\ = 9(?yi/5, with gA = 1-26 the axial coupling constant measured in neutron /3-decay. 
However, we will not use this relation but rather consider the value determined from the fit to the 
data as a check of the large N c expansion of QCD. In addition, there is the dimension two LEC 
ai. It only leads to a mass shift of the delta and can thus be absorbed completely in the physical 
nucleon-delta mass splitting. All these couplings are accompanied by separate off-shell parameters, 
which we (collectively) have denoted by zq, z and z'.# 6 These off-shell parameters are not observable, 
so they can be chosen freely. According to the common practice, we set z — z' — zq — |, which is of 
course a little bit arbitrary. As stated before, due to the explicit appearance of the A in the theory, 
there are new 1/m corrections to £ n N, which can be absorbed into the LECs q, di of the nucleon 
sector. This also means that a certain part of these LECs is explained by A properties (as it is well 
known, see Ref. [5]). In the following, we will not absorb these pieces in the q and di LECs (with one 

#6 Note that the special treatment of the leading nNA coupling is done for historical reasons. 
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exception to be given below), because such a redefinition also depends on the off-shell parameters 
z, z' and z . If one were to perform this redefinition, it would take the form 

c 2 — c 2 -2^(3 + 4. 2 )+8a|^(1 + 2. 2 ), (3.14) 

c 3 — c 3 + %^(1 + 8^ + 12^ 2 )-2A^(1 + 6z + 8z 2 ) , (3.15) 
9m 9m z 

c 4 — c 4 + %^(1 + 8^ + 12z 2 )-2A^f(l + 6^o + 8^ 2 ) , (3.16) 
(di + d 2 ) — (di + 4) + i ^- r (5 + 42 ;o )-2 — Zo --___(l + 4z + 12^o) , 



(3.17) 

), 

(3.18) 
(3.19) 

(rfi4-rfi 5 ) — ► {d 14 -d 15 ) - — -(l + 4z ) -8 — z -2—q—(1 + 4:Z + 12zzo). 

(3.20) 

In this case, the q now also have contributions proportional to A/m. Since c 2 and c 4 also appear 
in the third order amplitude, using these renormalized LECs would lead to contributions of fourth 
order, which is beyond the accuracy of our calculation. This we consider another reason not to make 
this redefinition. Similar problems related to the mixing of various chiral orders arise also in the 
fourth order analysis performed in the pure pion-nucleon EFT, see Ref.[2]. However, in one case, 
we have to perform this redefinition. The contribution proportional to b$ has to be absorbed into 
(di + d?), ds, g?5 and (du — J15). The resulting LECs and combinations thereof are labeled by a 
subscript "A". If that is not done, one introduces an additional redundant fit parameter. So this 
is the only redefinition we are making. The other effects from 6g can then be absorbed into the 
following combinations 

(/i + / 2 )-^ and (2/ 4 -/ 5 )-^-. (3.21) 

So we end up with the nine LECs from C^ffl plus four LECs (or combinations thereof) from C*%n\ 
(counting the leading 7riVA coupling as a free parameter, although we also perform fits with its value 
fixed) and one LEC from the An sector. 

4 Results 

4.1 The fitting procedure 

There are various possibilities to fix the LECs. We proceed here along the similar lines as in Refs. [1,2], 
namely we fit to the phase shifts given by different partial wave analyses in the low energy region. 
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This allows for a better comparison with the results obtained in the chiral expansion. As input we use 
the phase shifts of the Karlsruhe (KA85) group [14] and from the analysis of Matsinos [15] (EM98). 
Of course, there is also the phase shift analysis of the VPI/GW group, which we do not use here for 
two reasons. First, the solution called SP98 from the VPI/GW group [16] is no longer available (it 
was used in Refs. [1, 2]) and, second, no threshold parameters for the newest solution SPOO have been 
published. Clearly, once these are available we can update our investigation, but we do not believe 
that it will lead to wildly different results and new insight compared to what is presented below. 
Since no uncertainties are available for the KA85 solution, we mimic the ones of the Matsinos analysis 
in that case, which is 1.5% for S 3 i, 0.5% for Su, 1% for P 33 and 3.5% for the other P-waves. This 
assignment gives more weight to the better determined larger partial waves and is more natural than 
one common global error. We remark that the Matsinos analysis only includes data up to 200 MeV 
pion laboratory momentum. For this analysis, we proceed as in Refs. [1, 2], namely fit to the S- 
and P-wave phase shifts for momenta between 40 and 100 MeV and predict the phases at lower and 
higher energies. For the KA85 solution, we extend the fit region up to 200 MeV and thus predict the 
phases up to about 300 MeV (and of course also in the threshold region). As in Refs. [1, 2] the LEC 
dis is fixed by means of the Goldberger-Treiman discrepancy, i.e. by the value for the pion-nucleon 
coupling constant extracted in the various analyses. The actual values of g n N are g n ^ = 13.4 ±0.1 
and 13.18 ± 0.12 for KA85 and EM98, respectively. For the terms including the A, we perform fits 
with g w NA fixed (at the value of 1.05) and letting it free (as discussed above). Throughout, we use 
g A = 1.26, F n = 92.4 MeV, m = 938.27 MeV, A = 294 MeV and M w = 139.57 MeV. 

4.2 Phase shifts and threshold parameters 



LEC 


Fit 1 


Fit 1* 


Fit 2 


Fit 2* 


Cl 

c 2 
c 3 
c 4 


0.77 ±0.02 
-17.9 ±0.03 

20.2 ±0.06 
-15.6 ±0.04 


0.39 ±0.02 
-15.5 ±0.03 

16.7 ±0.04 
-12.5 ±0.03 


-0.44 ±0.01 
-0.67 ±0.03 
-0.07 ±0.03 
-2.51 ±0.04 


-0.32 ±0.01 
-1.59 ±0.03 
1.15 ±0.03 
-3.44 ±0.04 


(di + d 2 )A 
(4)a 

(4)a 

(g?14 — G?i 5 )a 


-5.91 ±0.06 
7.68 ±0.06 
-1.07 ±0.03 
-5.18 ±0.20 
-0.98 ±0.19 


-5.81 ±0.07 
6.60 ±0.07 
-0.59 ±0.04 
-0.09 ±0.21 
-0.97 ±0.19 


0.03 ±0.04 
-0.93 ±0.04 

0.75 ±0.02 
-0.44 ±0.15 
-1.35 ±0.14 


-0.36 ±0.04 
-0.44 ±0.04 
0.71 ±0.02 
-0.60 ±0.15 
-1.38 ±0.14 


g-xNA 

h + h 

f _|_ f Al 
J 1 ~ J 2 2m 

2/4 — fb — 4^ 
91 


1.32 ±0.03 
-12.0 ±0.10 
29.3 ±0.30 
40.2 ±0.73 
-1.42 ±0.02 


1.05* 
-11.1 ±0.13 

32.4 ± 1.30 

53.5 ±0.92 
-2.65 ±0.03 


0.98 ±0.05 
0.51 ±0.09 
-19.1 ±0.57 
-30.3 ±0.92 
-1.10 ±0.04 


1.05* 
-0.28 ±0.08 
-18.6 ±0.51 
-27.0 ±0.82 
-0.94 ±0.04 



Table 1: Values of the LECs in appropriate units of inverse GeV for the various fits described in the 
text. 
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We are now in the position to present results. For the analysis of Matsinos, we use 17 points for 
each partial wave in the range of q w = 41.4 — 96.3 MeV. We have performed fits leaving the coupling 
9itNA free or fixing its value at 1.05. We call these fits 1 and 1*, respectively. Therefore, we have 
to determine 14 and 13 parameters for these two fits. Since the data basis of Ref. [15] includes only 
data up to 200 MeV, we use these fits mostly for comparison with the results of the chiral expansion 
presented in Refs.[l, 2]. This is different for the Karlsruhe data basis, which extends up to very high 
energies. Consequently, for this case (KA85), we have fitted to the data up to 200 MeV pion lab 
momentum (i.e. 10 points per partial wave at q n = 40, 60, 79, 97, 112, 130, 153, 172, 185, 200 MeV). 
We extend the fits to higher energies than it was done in Refs.[l, 2]; this is possible due to the explicit 
inclusion of delta and allows to study explicitly the resonance region. Again, we perform fits with 
varying and fixed g n NA, which are denoted by fits 2 and 2*, respectively. The resulting LECs are 
collected in table 1. We point out again that the values for the LECs from the nucleon sector (the Cj 
and di) can not be compared with the ones obtained in the chiral expansion for the various reasons 
discussed above. For the Matsinos data, the x 2 /dof is very small and slightly better for fit 1 than for 
fit 1*. The reasons for this very small x 2 /dof are discussed in some detail in Ref. [2]. For the KA85 
phases, the x 2 /dof is essentially the same for both cases (and larger than for fits 1,1*, see again 
Ref. [2].). We also note that while the value of g n NA is somewhat larger than the canonical value of 
1.05 for fit 1, it comes out slightly lower in the case of fit 2. In the following, we will mostly discuss 
the results of the fits 1*,2* with g^NA — 1-05. Most LECs come out of natural size, i.e. of order one, 
with the exception of the LEC combinations from 0(e 3 ). This can be traced back to the fact that 
there are large cancelations at this order between the loop graphs with intermediate deltas and the 
counterterms. Such a phenomenon was also observed in the third order SSE calculation of neutral 
pion photoproduction, see [17]. Also, in case of fits 1,1*, the values for the dimension two LECs are 
quite large, and, in particular, the value for c\ is positive. Furthermore, it is important to stress 
that we find some sizeable correlations between the LECs C\^,z- This is not unexpected since the 
corresponding terms contribute only to the small isoscalar S-wave scattering amplitudes (in certain 
combinations, e.g. the S-wave isoscalar scattering length is only sensitive to the LEC combination 
— 2ci + C2 + C3). This will be taken up in the next subsection. We note that the value for the 71 A 
coupling gi comes out very different from the large N c prediction gi = 9g^/5 = 2.27. Note, however, 
since g\ only appears in the third order loop contribution, one can not expect to pin it down very 
precisely. 

The resulting fits and predictions based on the EM98 and KA85 phases are shown in Figs. 3 and 
4, respectively, in comparison to the results based on the third and fourth order chiral amplitudes. 
The number of LECs to be fitted was 9 and 14 in these cases. The results of the third order SSE 
calculation are clearly better than the ones based on the third order chiral expansion and comparable 
to the fourth order results, although the overall description is still slightly better in the latter case. 
As expected, the most prominent improvement can be found in the P33 partial wave, which is well 
described up to the delta pole, see Fig. 4. At lower energies, one can predict the threshold parameters, 
which are collected in table 2 for the fits 1* and 2* in comparison with the direct determination from 
the Matsinos and Karlsruhe phase shift analyses. The overall consistency is satisfactory but not as 
good as in the case of the fourth order chiral expansion. Again, we find a significant improvement of 
the scattering volume in the P33 channel, as expected due to the explicit inclusion of the delta. 

It is also important to discuss the issue of convergence. For that, we redo the fits for the amplitudes 
at first and second order in the small scale expansion. To leading order, one only has the coupling 
SVjva, whereas at second order one has the additional four LECs from the nucleon sector and one LEC 
combination from C^ A . The resulting phase shifts are shown in Fig. 5 for the EM98 analysis and in 
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Obs. 


Fit 1 


Fit 2 


EM98 


KA85 


a o+ 


0.41 


-0.94 


0.41 ±0.09 


-0.83 


°o+ 


-4.22 


-4.60 


-4.46 


-4.40 


a 0+ 


7.74 


8.95 


7.73 ± 0.06 


9.17 


b o+ 


1.42 


1.63 


1.56 


0.77 


af_ 


-5.49 


-5.54 


-5.46 ±0.10 


-5.53 


af + 


13.08 


13.27 


13.13 ±0.13 


13.27 




-1.21 


-1.46 


-1.19 ±0.08 


-1.13 




-8.21 


-8.14 


-8.22 ±0.07 


-8.13 



Table 2: Values of the S- and P-wave threshold parameters for the fits 1* and 2* in comparison 
to the respective data. The results for fits 1 and 2 are similar and thus are not given. Units are 
appropriate inverse powers of the pion mass times 10~ 2 . 



Fig.6 for the KA85 case. While the first order result is only good in the P 33 partial wave, the second 
order fits are of comparable quality than the third order ones, although at third order the x 2 /dof 
is significantly better. Such a behavior does not come completely unexpected since one expects the 
delta Born graphs to play the most significant role. Such an expectation is build on experience 
with many models that include the delta or also the explicit calculations of Compton scattering off 
nucleons in the SSE [18]. This behavior is different from what is found in the chiral expansion, where 
one still has large corrections when going from second to third order but mostly modest ones from 
third to fourth order, see [1, 2]. This means that in the channels where the delta plays a significant 
role, the resummation of higher order terms in the chiral expansion is important and well described 
in the approach used here. It is also interesting to study the convergence of the S-wave scattering 
lengths, as it has been done for the chiral expansion in [2]. In the small scale expansion we obtain 
the results collected in table 3. The convergence for the isovector scattering length is similar to what 
is obtained in the chiral expansion (as expected from the arguments presented in Ref. [19]). The 
isoscalar S-wave scattering length receives a large correction when going from second to third order, 
indicating that certain fourth order pieces related to pion-nucleon physics are still missing (as can 
also be inferred from the study in Ref. [2]). Some of these results were also found by Ellis and Tang [7], 
although their approach is based on a relativistic treatment of the fermion fields. Therefore, in their 
approach all 1/m corrections are resummed, but the delta is treated in a less systematic manner. 
Datta and Pakvasa [20] had already used the one loop representation of Ref. [21] and also added the 
delta, but not in a systematic fashion as done here. They also found a much improved description 
of the P 33 partial wave. 

4.3 The sigma term 

So far, we have exclusively considered the amplitudes in the physical region. We have noted that 
there are some strong correlations between some LECs. To further address that problem, we need 
additional input. This is provided e.g. by the so-called pion-nucleon sigma term, which is the 
expectation value of the QCD symmetry breaking terms in a proton (or neutron) state. It can be 
derived from the scalar form factor of the nucleon, 

<j(t) = {N(p')\m(uu + dd)\N(p)) , t = (p'-p) 2 , (4.1) 
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U(e) 


U(£ ) 


U(e ) 


4 + EM98 


0.0 


0.24 


0.41 


KA85 


0.0 


0.44 


-0.94 


a 0+ EM98 


7.90 


7.90 


7.74 


KA85 


7.90 


7.90 


8.95 



Table 3: Convergence of the S-wave scattering lengths. 0(e n ) means that all terms up-to-and- 
including order n are given. Units are 10~ 2 /M n . 



with \N(p)) a nucleon state of four-momentum p and rh the average light quark mass. The sigma 
term is nothing but the scalar form factor at t — 0, a = a(t = 0). We remark that the sigma term 
can not be obtained directly from scattering data. One usually considers its value at the Cheng- 
Dashen point, t = 2M 2 , where the chiral corrections are minimized. E = a(t — 2M 2 ) and a differ 
by 17 MeV, with 15 MeV stemming from the scalar form factor [22] and (at most) 2 MeV from the 
so-called remainder [23]. In addition, there exists a whole family of relations between E and certain 
combinations of threshold parameters, as detailed in Ref . [24] . These relations have been worked out 
to third order in the chiral expansion. We will use here the version given in Ref. [25] , 



E = ttF 2 [(4 + 2/i + /!>+, - AMX+ + 12/iM>+ + ] + E , 



(4.2) 



with E = —12.6 MeV and p = M n /m ~ 1/7. A special variant, which also contains some fourth 
order pieces, has recently been given by Olsson 



E = [F*F(2M 2 n )}, 
F(2Ml) = 14.5 a+ + - 5.06 {a^f - 10.13 (a^ 2 ) 2 



(4.3) 



16.65 6J + 



0.06 a+_ + 5.70 a+ + -0.05 , 



with the quantities on the right-hand-side being given in units of the pion mass. We will use these 
sum rules to further constrain our LECs. 

To third order in the small scale expansion, the scalar from factor of the nucleon reads (the scalar 
sector has also been discussed by Kambor [27]): 



a(t) 



+ 



4M 2 Cl 

9nNA 



m 2 


3M 2 




8 


2M 2 
3 


2J (- 



-2J (0) + (t-2M 2 )K (t,0) 



2 J (-A) + (2M 2 - t - 2A 2 )K (t, -A) - 2AI (t) + 



A 



8tt 2 



(4.4) 



in terms of the standard loop functions listed in [21]. For t — 0, this gives 
-4M 2 Cl - 



a(t = 0) 



g A \ 2 9M 3 fg nNA 



64tt 



2 m 2 , / A JE 2 



M 2 



M 



(4.5) 



Note that the terms ~ g^NA m Eqs. (4.4,4.5) are of fourth order in the chiral expansion, as already 
pointed out long time ago [28]. We note that a only depends on the LEC c\ (and also the coupling 



17 



9nNA f° r the where it is left free). If we now use C\ as determined in fits 1 and 1*, we get an 
unphysical negative sigma term because c\ is positive and large. That, however, is an artifact since 
not all the q's can be determined independently. For the fits 2 (2*), we obtain the following values, 



a = 51.1 (47.3) MeV 



(4.6) 



These numbers are slightly larger but consistent within error bars with what has been found before, 
see Refs.[22, 29]. However, one might question the accuracy of this determination because also in 
this case one has large correlations between certain LECs. To overcome this, we have performed a 
different set of fits. As additional input we take the sigma term as determined from the sum rules, 
Eq.(4.2) and Eq.(4.3) using the threshold parameters from the EM98 and the KA85 analysis as input. 
More precisely, since the two sum rules differ by some terms of fourth order (and higher) in the chiral 
expansion, for each set of input threshold parameters we get two numbers for a = S — 17 MeV, which 
we average to obtain the central value and their spread is taken as the theoretical uncertainty. This 
gives a = (58.5 ± 5.4) MeV for EM98 and a = (45.5 ± 2.7) MeV for KA85. The corresponding fits 
with this additional input are denoted by fit 1* and fit 2*. The resulting LECs are shown in table 4. 



LEC 



Cl 

c 2 
c 3 
c 4 



(d l + g? 2 )a 
(4)a 
(4)a 
{du — di 5 )A 
<ii8 



9nNA 

h + b 8 

fl + J 2 ~ ^ 

2/ 4 - h - 
9i 



Am 



Fit 1' 



-0.18 ±0.02 
-5.72 ±0.03 
6.05 ±0.03 
-8.93 ±0.04 



5.52 ±0.07 
-4.12 ±0.07 
-0.89 ±0.04 
-18.8 ±0.22 
-0.99 ±0.20 



1.27 ±0.04 
-7.33 ±0.12 
-31.3 ±0.98 
-68.0 ±0.77 
-2.05 ±0.02 



Fit 2* 



-0.35 ±0.09 
-1.49 ±0.66 
0.93 ±0.87 
-3.08 ±0.81 



-0.57 ±0.64 
-0.30 ±0.64 
0.74 ±0.10 
-0.10 ± 1.60 
-1.34 ±0.24 



1.00 ±0.08 
-0.03 ±0.89 
-17.9 ±4.49 
-23.0 ±9.78 
-1.05 ±0.41 



Table 4: Values of the LECs in appropriate units of inverse GeV for the fits using as additional input 
the sigma term. 



First, as expected from the numbers given in Eq.(4.6), there are only minor changes in case of fit 2* 
as compared to fits 2,2*. This is very different for fit 1* compared to fits 1,1*. The value of C\ is 
now negative and also, the pion-nucleon dimension two LECs have more natural values. In addition, 
the correlations between the LECs 0^2,3 are somewhat smaller than before. A very important result 
is the stability of the coupling constant g n NA, which comes out consistent with what was found in 
fits 1 and 2, respectively. 
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A Pion-nucleon-delta vertices 



Here, we give the Feynman rules for the relevant 7riVA vertices. Consider first the case on an 
incoming nucleon (momentum pi), an outgoing delta (momentum P M , isospin i), and an outgoing 
pion (momentum q, isospin a): 
1st order: 

^V<r . (a.i) 



2nd order: 

3rd order: 

~„ 1 



^3 + ^8 ucia 9ttNA d ria 

-v ■ qq^o —P^v - go . 



mF 



(A.2) 



-Qn({ei + e 2 )(v ■ gf + 2M 2 (-2e 4 + e 5 



m 
m 



-l(v ■qv-(P + p 1 )-q-(P + pi)) - b 8 q ■ Pi - 



z{{vqf-q 2 ) + {h + b,){v.qf 



v • Pv ■ pi — P ■ pi) 



_9^fA ( _ 1 (u pv q _ p q) _ {2zq _ 1)v . Pv . q + {2zq _ A)Av . o 1 (A3) 
6m 2 J J 

Similarly, for the case of an outgoing nucleon (momentum p 2 ), an incoming delta (momentum P M , 
isospin i), and an outgoing pion (momentum q, isospin a), we have: 
1st order: 



2nd order: 



3rd order: 



+ b ~^v ■ qg »5» - 9 -^P,v ■ q5* 
F mF 



(A.4) 
(A.5) 



5 ia - 
F 



-<7 M ((ei + e 2 )(v ■ qf + 2M 2 (-2e 4 + e 5 )) 



■ qv ■ {P + p 2 ) - q - (P + P2)) + hq ■ p 2 



+ - 







m 


2 v 


Pa 


r h 


m 


L 2 ; 




flViVA 


3m 



9nNA 
Am 



v ■ Pv ■ p 2 - P ■ p 2 ) 



z((vqY-q 2 )-(b 3 + b s )(vqy 



'--(v -Pvq-P-q)- (2z - l)v ■ Pv ■ q + (2z - A)Av ■ q) 



(A.6) 
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Note that possible tvtiNA vertices only start at second order and can therefore not appear in the 
leading loop graphs of order 0(e 3 ). This explains the vanishing of some diagrams not drawn in Fig. 2. 



B Expressions for the 1/m corrections 

In this appendix, we give the lengthy expressions for the various 1/m corrections stemming from the 
A-insertions, which are of relevance to our problem. 

1/m corrections to C^n'- 



N lo B%hoC A m B$ A N = - 9 -f^N { 1(1 + 8^ + 12^)5.^5 

1 
3 



+ Uz-Zzq-AzDv-w^^vw^N . (B.l) 



1/m corrections to £^n : 



N { [7o(4 1 ^ 1(0) ^i)-So4 0) " 1 ^ + h.c.]+ 7 oBW 7o Ci 1) - 1 ^ 

+ [ 7 o£lS7oCi OM ^i + h.c.]}N 
a 2 2 

= N { f^£-[-(l + 4z + 12£)S ■w i ? j vw j 2iS-D + h.c. 
1 (2m) 2 L 3 

--(3 + 4z + Az 2 )v ■ w^S ■ w j 2iS ■ D + h.c. - -(1 + 4z + 4z%)S ■ w'fHv ■ DS ■ w> 

3 3 

16 8 
+—(1 + 6z + 8z$)AS ■ w'f'S ■ w j + -(5 + 8z + 12z 2 )S ■ w^HS ■ Dv ■ w j + h.c. 

3 9 

2 .. . 

+ -(-1 + 4z - 4zq)v ■ w l iv ■ DC 3 v ■ w J + 4(1 - 2z )Av ■ w l C J v ■ u 

\l + Az + Az + 12zz )(b 3 + b 8 )S ■ w^K^SaVp - h.c. 



QnNA r [ 

lr 
4 



2m 



3 



+-(1 + Az + 12zz )b 3 S ■ w'C'ivipVaSp - h.c. 
-|(6 + 13z + Azz )b 3 v ■ w^KpSaSp - h.c. 



+ I(_5 + 4 ^ + 4^ + Azz )(b 3 + b 8 )v ■ w^w^vp - h.c] }N . (B.2) 
1/m corrections to C^na- 



1/m corrections to C^^. 



T ul0 B A {1) C A m B$ A N + h.c. = - 9 -^2f i u iD i ^ kl v ■ w l N + h.c. . (B.3) 
^ 2m 

(3) 
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TlboVTh& m B$ + ToB* 1 + ToB* 1 W <0) B« }jV + h.c 

■ h ' A ^iS • D ik £ kj wfciS ■ D - 8 + 32z ° iS ■ D ik £ kl iD l ™£ mn S ■ w r 



3 



(2m] 

+ ^z iD ik £ kl iS ■ D lm C m S ■ w n - 8Z °~ A iD lk i kl iv ■ D lm C n v ■ w n 
+ 8Zo J 16 AiD^v ■ w 1 } }N + h.c. . (B.4) 



3 

1/m corrections to C nA . 

1 



T^oBf^oC^B^Ti = —T^^S-Dg^S-Tl. (B.5) 



1/m corrections to C^ A : 

T^oBf^oC^B^Tl = -j^Tftv {2iS ■ Dg^iv ■ D + A)2iS ■ D 

+4iD^iv ■ D - A)iD u }Tl . (B.6) 

C Tree and loop amplitudes 

In this appendix, we give the lengthy analytical expressions for the tree and counterterm as well as 
one loop amplitudes involving intermediate delta states, as depicted in Fig.l (tree and counterterm 
graphs) and Fig. 2 (loop graphs). We use the loop functions defined in Ref.[21]. All other notation 
has been defined previously. 

Tree diagram contributions: 

= -gl NA 2 -(2^-2M* + t) / ' ' 



'9 V ' * \uj - A u + A 

+ 9 l NA ±-{(2^ - 2MI + ^ - Ml) ((-l^) 2 + (-1 

-(2a. 2 - 2Ml + t)(4c 2 - 4M 2 + t) (^) 2 

+Auo{Auo 2 -4Ml + t)—!— 
v n 'uo + A 

-4cu 2 (3 + Az 2 ) + (2M 2 - 0(1 + 8^o + 12^ 2 )} 
- g nN A(b 3 + b 8 )^u(2co 2 -2Ml + t)(^— + 



2^ 

.1' 



9 v n J \u-A u + A 
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u + A 



+ l -(2u: 2 - 2M 2 + t) 2 (AuJ 2 - 4M 2 + (^^) 3 
-uj(Auj 2 - AM 2 + t)(8uj 2 - 8M 2 + 3f) (j^) 
+ I(4cu 2 - 4M 2 + t)(72cu 2 - 24M 2 + 

+ 16Acj 2 (1 + 2zl) - 2A(2M 2 - t){\ + 6z + 8^) - 2uo{Auo 2 - AMI + *)( 3 + Az 2 )} 
+ 9«NA(h + b 8 )^-{2uj(uj 2 - M 2 )(2uj 2 - 2M 2 + t) 

+2u(2u 2 - 2M 2 + t)(Auo 2 - AMI + *) {j^£ 

-2[{u 2 - M 2 )(22cu 2 - 6M 2 + 5f) + t 2 ]—^— 

u> + A 

-2(^ 2 - M 2 )(2^ 2 - 2M 2 + Oj^} 
- ^A&8^-(2^ 2 -2M 2 + t)(cu 2 -M 2 ) 



uj-AJ W + A 



9m 1 " /v Vcj- A uj + A/ 
+ [2^ A (ei + e 2 ) - (6 3 + hf]\uj 2 {2uj 2 - 2M 2 + t) (—^—r ^ 

- ^ A (2e 4 - e 5 )^M 2 (2^ 2 - 2M 2 + t) - . (C.l) 

o 2/1 1 \ 



+ (4w 2 -4M 2 + t) (— 



A/ Vcj + A, 
2 



-4u,^-} 
cj + A- 1 

+ 57tJVa(&3 + h)-v ( r- — r) 

9 \c<j — A uj + A/ 



+i(2^ 2 - 2M 2 + t)(4^ 2 - 4M 2 + t) (^) 3 
-lo(10lo 2 - 10M 2 + St) ( —1 -£ j 

+ (16u 2 - 12M 2 + 3t)—^— 
v n j uj + A 



22 



-2U(1 + A Z 

+ 9nNA(h + b 8 )^-{ -2uj{uj 2 -MI) 

+ 2^-4M^ + t)(-l^) 2 
-2(7^-3M 2 w + t)—L- 

UJ -\~ L\ 



uj-AJ \lu + A 



1 



— 16uz 

+ 9^h^{{u 2 - Ml) + _L_) _ Ml + 4, + 12**)} 



[2^A(ei + e 2 )-(6 3 + 6 8 ) 2 ]-o; i 



2 2 / 1 1 



-A + A 



u-AJ W + A 

2 



+ ^A(2e 4 -e 5 )^(-^ + -^) . (C.2) 

-i(2^ 2 - 2M 2 + t)(4^ 2 - 4M 2 + (^)' 
+2^ 2 -4M 2 + t)-^} 

+ ^3+M^ 2 -2M 2 + t)(-^--^; 

+i(2^ 2 - 2M 2 + i)W - 4M 2 + t) 

- l -w(Au 2 - AMI + t)(8^ 2 ~ 8M 2 + 3t) (— ^ 

+ ^-(Alu 2 - AMI + t)(72iu 2 - 24M 2 + 7t)—^— 
L2 lu + A 

+w[2(2M 2 - 0(1 - 2^o - 2^ 2 ) - 12cu 2 + M 2 (9 + 4z + I2z 2 )}} 
+ ^a(6 3 + &s)^{ - w(2lu 2 - 2M 2 + t){oo 2 - Ml) ((^^) ? 

+w(2u 2 - 2Ml + t)(Auu 2 - 4M 2 + t) (—^-A 

\uj + A/ 

+ [(^ 2 - M 2 )(-22^ 2 + 6M 2 - 5t) - t 2 ]-^_ 



2 

+ 



+ A 
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+ {2u 2 -2Ml + t){u 2 -Ml)-^ 

-Auj{2Ml - t)zo + cj 3 (16 + 5z - Azz ) - luM 2 (Q + 13z + 4zz )} 
+ 9^h±{(^ - 2Ml + t)(w > - Ml) {-L- + _L^) 

-u{2Ml -t)(l + Az+ 12zz ) - cj 3 (4 + 5z - 2Qzz Q ) + w^ 2 (6 + Viz + Azzq)} 
[2^jvA(ei + e 2 ) - (6 3 + 6 8 ) 2 ]^ 2 (2^ 2 - 2M 2 + t) ' 

A / 1 



L *"'"" ' !) ^ ' \oo- A + u + AJ 

+ ^A(2e 4 -e 5 )^(2c, 2 -2M w 2 + t)(-^ + : ^) . (C.3) 



+ A, 



2 



+2u; 



uj + A 
+(1 + 8^0 + 12^)} 
2 / 

g^NAih + 6 8 )-w ( 4 



cj-A w + A, 



a; Vc^ + a, 

3 



+i(2w 2 - 2A/| + t)(4u, 2 - 4Af„ 2 + t) 

-r (10w2 - 10M - +3t) (jTA) 2 

+ l (18w » - 12M| + 3*)^ 
+u(l + 8^o + 12^q) - 2A(1 + 6z + 8^)} 

+ g*NA{h + bz)-^{u{u 2 -M\ 

+^-4M 2 +t)(-^; 



uo-A) Vcj + A, 



_ (7 ^_3M, 2 + i) — 

1 

'a; - A 

i\ I . 

\LU — A uo + A, 



uj — 

9nNAb s ^(u; 2 - Ml) 
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+ PWei + e 2 ) - (6, + h?] l r > - ^) 



Loop diagram contributions: 



= (JoH + JoH){ - 9? y^ 2 , (2M* - 4cc 2 M 2 - tM 2 + 2a; 4 + c 2 t) 



27(u; 2 - A 2 ) 



+ 243( 8 J 2 ^ A A2)2 (8MV + 10M 4 A 2 - 16M„V - 4M> 2 * 

-20M 2 u; 2 A 2 - 5M 2 tA 2 + 8cj 6 + 4uA + 10cu 4 A 2 + 5cu 2 tA 2 )} 

+ Jo(-A){ 2i87 ^J 2 A _ A (-972MV - 972M> 2 A 4 + 1944M> 4 A 2 

-3888cj 4 tA 2 + 1944cu 2 tA 4 + 1944cA) 

+ ^T^^^(" 792M > 2A2 - 504M> 4 + 7776M 2 c/A 2 
2187u; z (u; z: — A Z J Z 

+252M 2 cA + 396M 2 w 2 tA 2 - 5688M 2 cj 2 A 4 + 504M 2 u; 6 

-3492uAA 2 + 2844cu 2 tA 4 + 5688c/ A 4 - 6984uM 2 ) 

2 2 

+ 010 /: Af f 1 A2 , 2 (-50MV + 500M> 2 A 2 - 450M 4 A 4 + 50MV 
zlo7u; z (u; z — A^J Z 

-450M 2 cj 4 A 2 + 225M 2 tA 4 + 25M 2 uA - 250M 2 cAA 2 - 50M 2 cj 2 A 4 

+450M 2 A 6 + 500cu 4 A 4 - 225t A 6 + 250uA A 4 - 50uM 2 - 450cj 2 A 6 

-25uAA 2 ) 

50gl NA g Agi M 4 + 2 ^2 M2 + + _ 2 ^ 2a2 _ ^ 2 



243cj 2 

+ - ^AgL^ (_2MV + 18M 4 A 2 + 2M 2 cu 4 + M 2 tu 2 t - 1QM 2 lu 2 A 2 
27cj 2 (cj 2 - A 2 ) 

-18M 2 A 4 - 9M 2 tA 2 - 2c/A 2 - uHA 2 + 18cu 2 A 4 + 9tA 4 )} 

- J (A)—^^—(-2My - 2M 4 A 2 + M 2 tA 2 + 4M> 2 A 2 + 2M 2 A 4 
2([uj z — A A ) A 

+2M> 4 + M 2 oo 2 t - tA A - 2uo 4 A 2 - 2cu 2 A 4 - cu 2 tA 2 ) 
1 ( J {u-A) J (-u-A) \ ( 25g 2 nNA g 2 A2 . RM 4 A 2 A()M 2 ( a 
+ ^{(cu-Ay ( W + A) 2 ) i 4374 (2 ° M ^ +18M - A - 40M ^ 

-10M 2 cj 2 t - 96M 2 u; 2 A 2 - 9M 2 tA 2 - 18M 2 A 4 + 2(k> 6 + 78u; 4 A 2 
+10uH + 39cu 2 tA 2 + 18u; 2 A 4 + 9tA 4 ) 



259 * NA9A9l (4My + 2M 4 A 2 - 8M 2 c/ - 2M 2 u 2 t - 16M 2 u 2 A 



243 

-2M 2 A 4 - M 2 tA 2 + Ato 6 + Ucu 4 A 2 + 7u 2 tA 2 + 2cj 2 A 4 + tA 4 + 2uH) 



25 



+ %^(20M> 2 + 18M 4 A 2 - 40M> 4 - 10M 2 u 2 t - 96M 2 uj 2 A 2 
54 

-9M 2 tA 2 - 18M 2 A 4 + 20cu 6 + 78c/A 2 + 10w 4 t + 39w 2 tA 2 
+18cu 2 A 4 + 9tA 4 )} 
A ( Mu-A) J (-u-A) \ ( 25gl NA gf ^ ^ , 

+56M 2 A 2 + 10M 2 t - 60u; 4 - 56cj 2 A 2 - 30luH - 2UA 2 ) 

_ g5g|^4gl ( _ 4M 4 + 16M 2^2 + g M 2 A 2 + 2M 2 t _ ^4 _ 

-8cj 2 A 2 - 4tA 2 ) 

2 2 

+ ^ jV ^ A (-20M 4 + 80M 2 w 2 + 56M 2 A 2 + 10M 2 t - 60c/ 
54 



-56c/A 2 -30c/t-28tA 2 )} 



- J o(-A) 27 ^^ 2) A(-2M 4 + 2M 2 y + M 2 t + 2M 2 A 2 - 2c/A 2 - tA 2 ) 

- K (t, -A)^4p( 2M ' - 2M 2 A 2 - 5M 2 t + 2t 2 + 4tA 2 ) 

9 

+ J (t)^rAA(M 2 -2t) 



9 



+ 104976^ff-A^)^ (1296M - VA + 10368M ' A3 " 14 688M 2 A 5 



-5184M 2 tA 3 - 648M 2 oAA - 3456M 2 cu 2 A 3 - 1296M> 4 A + 14688a; 2 A 5 
-6912c/ A 3 - 3456c/tA 3 + 7344tA 5 ) 

fh ' A A(1944M 4 (u 2 - A 2 ) - 2916M 2 A 4 + 7776M 2 c/A 2 



104976tt 2 (c/ - A 2 ) 2 

+972M 2 tA 2 - 4860M 2 c/ - 972M 2 cA + 5832t(c/ - A 2 ) 2 ) 

+ 104976fff-A 2 ) 2 ( - 27 ° 07rM - ( " 2 + A2) " 3030M ' A ^ 2 " A2 ) 
+ 1350ttM 3 (cj 2 + A 2 ) (2c/ + t) + 1515M 2 c/tA - 9840M 2 A 5 
+ 18410M 2 c/A 3 - 1515M 2 tA 3 - 8570M 2 c/A - 21440c/ A 3 + 9840c/ A 5 
+4920tA 5 + 11600c/ A - 10720c/tA 3 + 5800c/tA) 
^91na9a9i 



:(20ttM*(i/ + A ) - 42M 4 A(c/ - A ) 



1944tt 2 (c/ - A 2 ) 2 

-10ttM 3 (c/ + A 2 )(2c/ + t) - 21M 2 tA 3 + 21M 2 c/tA - 112M 2 A 5 
-70M 2 c/A + 182M 2 c/A 3 + 112c/ A + 56tA 5 + 112c/ A 5 - 112c/tA 3 
-224c/A 3 + 56c/tA) 

+ 648 jg?! i A2)2 (-vrM 5 (42^ 2 - 150A 2 ) - 24M 4 A(c/ - A 2 ) 

+ttM 3 (21c/ - 75A 2 )(2c/ +t) + 200M 2 c/A 3 - 80M 2 c/A - 120M 2 A 5 
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+12AfVf A - 12M 2 tA 3 + 104c 6 A + 52tc 4 A - 112tc 2 A 3 + 60tA 5 

-224c 4 A 3 + 120c 2 A 5 ) . (C.5) 



+ 



+ j °(- A ){ 2187 !f ( i" A -\ 2 ) 2 (- 27oM v + 9i8 - 4 - 648 - 2a2 ) 

+^^¥^^( 50M '(^ 2 - A2 ) + l25A4 - 2 °°^ 2a2 + 75 ^ 4 ) 

zl87c(c z — /a z ) z 

- ^tsT 9 ^ + 27^ A -A^) (8M ' V " 4M ' A2 + ^ " 5A4) ) 
1 / J (c-A) J (-c-A) \ ( 25gl NA gf ~ 2 3 

+ wU^-A) 2+ (oT+Ap Ji 4374 A - 30c A - 32A ) 



25 £na9a9i {4M 2 A _ 6uJ 2 A _ 8A 3 } 

■NAi 

27 



243 

+ ^V^i (10M 2 A _ 15cJ 2 A _ 28A3) | 



1 / J (c - A) J (-c- A) \ ( 25gl NA gl , M 2 2 fi/w -2 A 2 , . ; 4 

+51c 2 A 2 + 6A 4 ) 
_25^Mi y _ 2a2 ^ 4 n ^ A2 2A4 
243 

+ ^g^i(-5M 2 c 2 - 18M 2 A 2 + 5c 4 + 63c 2 A 2 + 18A 4 )} 

+ J {-/\)iM n Aj| 72g ^ + 81(cJ 2_ A 2)- 243tc; 

, 2 91na9a \ 
3c J 

4 

_ 9ttNA IjJ A 2 C-2916M 2 4- 1944A 4N I 

26244vr 2 (c 2 _ A 2 ) 2 1 ^ ibM ^ + i9MA j 

- , f 2 M f a2 , 2 (1350ttM 3 A + 645M 2 (u> 2 - A 2 ) + 1700c 4 + 1860A 4 
524887r 2 (c 2 - A 2 ) 2 



-3560c 2 A 2 } 
91na9a9iuj _ M 2 



547r 2 (c 2 - A 2 )' 

,2 



+ 9 / 2 ' 2 fff A2)2 (-10^M 3 A + 3M 2 (c 2 - A 2 ) + 12(c 2 - A 2 ) 2 ) 

27 



+ 3247^1^ - A^ {nM ^ 8uj4 ~ 93u2A2 + 36A4) " 42M > 2A ^ 2 " A ') 



+4c/A 3 - 10u 2 A 5 + 6Au 6 ) . (C.6) 
F£g-(u,,t) = (JoH + J o(-^))^^{ gl( i y^ 2) ^A(-2M, 4 + 4M^ 2 + M 2 t - ^ - 2c 4 ) 

2^ata4(2M 4 - 4M> 2 - M 2 t + 2c/ + uH)} 



27 

+ (JoH - M-u))^^ 2 { 24 3^ A A2) ( 2 ^ A2 + 4 ^V - ^ 

-4M 2 cu 2 A 2 - 2M 2 w 2 t - M 2 tA 2 + 2c/ A 2 + 4cu 6 + 2cA + c/tA 2 ) 

_ 2 9lNA9A A (_ 2M 4 + 4M 2^2 + M 2 t _ _ ^ j 

+ J (A) ^f^^ 9 cc>A(-2M 4 + 2M> 2 + M 2 t + 2M 2 A 2 - 2c/A 2 - tA 2 ) 

+ Jo( " A H^||^i^ A (- 100M '^ 2 - A2 ) + 50M ^ 2t - 250M ' A4 

-50M 2 cj 4 - 50M 2 tA 2 + 300M 2 c; 2 A 2 + 125tA 4 - 400c/A 2 
+75uH + 250cj 2 A 4 - 200cu 2 tA 2 + 150cu 6 ) 

?tt7VA A /K/in A/T4, ,2 , ioncj/2, ,2 a 2 n^m\/r2, ,2 



2187^- A 2)2 A ( 540M V + 1296M> 2 A 2 - 270M 2 cu 2 t 

-2376M 2 c/ + 918c/t - 648c/tA 2 - 1296c/A 2 + 1836c/) 
^Ogl NA g A gi A/ nA , 2 , , .2 



- A(-2M 2 + 2c/ + t) 



- -"r4^(-16MV + 8M 4 A 2 + 16M 2 c/ + 8M 2 c/t + 10M 2 A 4 



243w 
^91na9a 
27wA(w 2 - A 2 ) 

-10M 2 c/A 2 - 4M 2 tA 2 + 2c/A 2 + c/tA 2 - 10c/ A 4 - 5tA 4 )} 
+ ~ ( MUJ ~ A ) + M ~ UJ :, A) ) ( 2% l^ l2 (-40M 4 + 100M 2 c, 2 + 64M 2 A 2 



uj 



(cj - A) 2 (uj + A) 2 ) I 8748 

+20M 2 t - 60c/ - 30c/t - 64c/A 2 - 32tA 2 ) 

_ 25^ jvA ^ 1 ( _ gM 4 + 2QM 2 u2 + AM 2 t + 16M 2 A 2 _ 12cJ 4 _ 16cJ 2 A 2 

486 

-6uH - 8tA 2 ) 

2 2 

+ 9jrNA9A ,_ 40M 4 + ioomW + 112M 2 A 2 + 20M 2 t - 60c/ - 112c/A 2 

j^Qg V 7T 7T 7T 7T 



-30c/t-56tA 2 )} 



-114M 2 c/A 2 - 5M 2 u 2 t - 6M 2 tA 2 - 12M 2 A 4 + 10c/ + 5uH 

.4 A 2 , r-|. .2, A 2 , 10, .2 A 4 , c , A 4\ 



+102c/A 2 + 51c/tA 2 + 12o/A 4 + 6tA 4 ) 



28 



259 " na9a9i -(2M^ 2 + 4M 4 A 2 - 4M 2 c 4 - 26M> 2 A 2 - M 2 w 2 t 



486 

-2M 2 tA 2 - 4M 2 A 4 + 2c 6 + uH + 22c 4 A 2 + 4c 2 A 4 + llc 2 tA 2 
+2tA 4 ) 

+ gIiVAgi (10M 4 2 36M 4 A 2 _ 20 M 2 c 4 - 5M 2 cu 2 t - 162M 2 c 2 A 2 

-18M 2 tA 2 - 36M 2 A 4 + 10c 6 + 126c 4 A 2 + 5c 4 t + 63c 2 tA 2 
+36c 2 A 4 + 18tA 4 )} 

j5(-a){ (^|^ + gff^ ) (a 2 - *w - 2- 2 - o 

9 

509l NA 9A9l ( _ 2M 4 + 2M 2^2 + M 2 t + 2M 2 A 2 _ ^2^2 _ ^2) 



243c 

+ ^VA^ ( _ 2M 4 + 2u 2 M 2 + M 2 t + 2M 2 A 2 _ ^2^2 _ ^2)1 

3c ^ 

tf (t, -A)^4r^ A ^(- 2M ' + * + 2A 2 ) 

9 

h{t)^4^uo{-ml + 12A 2 + 5t) 

209952^p-A^)^ (1944M - (a;2 " A2) " 5832M ^ 2A2 " 972M 2 ^ 2 t 
+ 1944M 2 c 4 + 972M 2 tA 2 + 3888M 2 A 4 - 2268t(c 2 - A 2 ) 2 ) 

2 2 

;/:; x xm -c(-54007rM*A - 3030M 4 (c 2 - A 2 ) 



209952tt 2 (c 2 - A 2 ) 2 

+27007rM 3 A(t + 2c 2 ) - 3790M 2 cj 2 A 2 + 1260M 2 A 4 + 2530M 2 c 4 
-1515M 2 tA 2 + 1515M 2 c 2 t + 500c 6 + 760c 4 A 2 - 1260c 2 A 4 
-630tA 4 + 380c 2 tA 2 + 250c 4 t) 

7 c(16200M 4 A 2 - 4536M 4 c 2 + 22Q8M 2 uj 2 t 



209952tt 2 (c 2 - A 2 ) 2 v 

-5400M 2 cj 2 A 2 - 18576M 2 A 4 - 8100M 2 tA 2 + 4536M 2 c 4 

.2. a 2 , nr,oojA4 i nonn, A A 2 , i orr^c. .2 A 4\ 



-5400c 2 tA" + 9288tA 4 - 10800c 4 A" + 18576c 2 A 4 
+ 3 88 3^ff A2)2 ^(40 7 rM|A - 42M> 2 - A 2 ) - 20vrM^A(t + 2c 2 ) 
+M 2 (c 2 - A 2 )(21t + 14c 2 - 14A 2 ) + 14(c 2 - A 2 ) 2 (2c 2 + t)) 

1 oofi 2 l lf\ A2^ 2 ( 7rM '( 144A4 + 192w * ~ 372c 2 A 2 ) - 96M 4 c 2 A(c 2 - A 2 ) 
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-388cj 2 A 3 + 178cj 4 A) . (C.8) 

D Threshold parameters 

Here, we collect the analytical expressions of the A-contribution for the S- and P-wave scattering 
lengths and effective ranges. The corresponding pure nucleon terms are given in Ref. [2]. As before, we 
explicitly display the off-shell parameters but we note that these are fixed in the actual calculation. 
We have: 
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We note that some of these terms are clearly of fourth order in the chiral expansion, because in that 
case one counts A as order O(p ). 
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Figures 



Figure 1: Tree and counterterm graphs involving the delta. Dashed, solid and double lines refer to 
pions, nucleons and deltas, in order. Crossed graphs and diagrams that vanish are not shown. The 
vertex dot and vertex square refer to an insertion from the dimension two, respectively three effective 
7rA or nNA Lagrangian. 
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Figure 2: Leading one loop graphs involving the delta. Dashed, solid and double lines refer to pions, 
nucleons and deltas, in order. Crossed graphs and diagrams that vanish are not shown. 
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Figure 3: Fits and predictions for the EM98 phase shifts as a function of the pion laboratory mo- 
mentum q n . Fitted in each partial wave are the data between 41 and 97 MeV (filled circles). For 
higher and lower energies, the phases are predicted as shown by the solid lines. Dotted and dashed 
lines: Third and fourth order calculation based on the chiral expansion [1, 2]. 
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Figure 4: Fits and predictions for the KA85 phase shifts as a function of the pion laboratory mo- 
mentum q n . Fitted in each partial wave are the data between 40 and 200 MeV (filled circles). For 
higher and lower energies, the phases are predicted as shown by the solid lines. Dotted and dashed 
lines: Third and fourth order calculation based on the chiral expansion [1, 2]. 
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Figure 5: Fits and predictions for the EM98 phase shifts as a function of the pion laboratory mo- 
mentum q n to first (dotted lines), second (dashed lines) and third (solid lines) order in the small 
scale expansion. Fitted in each partial wave are the data between 41 and 97 MeV (filled circles). For 
higher and lower energies, the phases are predicted. 
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Figure 6: Fits and predictions for the KA85 phase shifts as a function of the pion laboratory momen- 
tum q n to first (dotted lines), second (dashed lines) and third (solid lines) order in the small scale 
expansion. Fitted in each partial wave are the data between 40 and 200 MeV (filled circles). For 
higher and lower energies, the phases are predicted. 
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